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A new structural model of Au induced Si(553) surface is proposed. The model accounts for recently 
experimentally found value of the Au coverage, i.e. 0.48 monolayer, which suggests formation of two 
gold chains on each Si(553) terrace. The resulting structural model, like the models of other vicinal 
Si surfaces, features the honey-comb chain, but there is no buckling at step edge, which is observed on 
Si(335)-Au and Si(557)-Au surfaces. The present model is more stable than the models with single 
Au chain only, and agrees very well with existing experimental data. In particular, calculated band 
structure, featuring two metallic bands coming from hybridization of the gold in both chains with 
neighboring Si atoms, perfectly matches the photoemission data. Moreover, theoretical scanning 
tunneling microscopy topography remains in excellent agreement with the experiment. 

PACS numbers: 73.20.At, 71.15.Mb, 79.60.Jv, 68.35.B-, 68.47.Fg 



I. INTRODUCTION 



More than decade ago it has been discovered that sub- 
monolayer coverage of gold stabilizes stepped silicon sur- 
faces, leading to one-dimensional ordering of the surface 
Since then various stepped Si surfaces have been ex- 
tensively studied in relation to expected exotic phenom- 
ena, characteristic for systems of reduced dimensionality 
[2, 3]. Perhaps one of the most known examples is a gold 
decorated Si(553) surface. The Si(553)-Au surface con- 
sists of Si(lll) terraces 4| x a^n^ = 1.48 nm wide and 
single atomic steps. The Si(553) surface normal is tilted 
from [111] direction towards [112] by 12.5° [4]. The prop- 
erties of this surface have been investigated by number 
of techniques, including scanning tunneling microscopy 

(stm) H H, @, 0, i, IE [TM3L a ^i4 e resolved P h °- 

toemissionfARPES) 0, S, 0, 0, El Q, x-ray diffrac- 
tionJS, [H ElL and density functional theory (DFT) 



Topography of the surface, as measured by STM, fea- 
tures single few nanometers long chains on each terrace, 
which have been identified as originating from the step 
edge Si atoms. Those chains are observed to have 3.84 
Aperiodicity along the chains, i.e. in [110] direction, 
which perfectly matches the Si lattice constant in this 
direction. This is completely different from what is ob- 
served on other vicinal Si surfaces, like Si(335)-Au 0, 21] 
or Si(557)-Au [3, [H, |23[ , where chains show periodicity 
doubling, i.e. STM topography shows maxima along the 
chains with a period of 2 x a^iop The periodicity dou- 
bling on Si(557)-Au and Si(553)-Au surfaces has been 
explained in terms of the buckling of the step edge Si 
atoms [24, 25] , according to which every second step edge 
Si atom occupies up (down) position. On Si(553)-Au sur- 
face one occasionally observes periodicity doubling and 
even oscillations of topography with a period of 3 x a^jo] 
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[8[ , but this has its origin in defects which force the chains 
to multiply their periodicity ^ n general, the ob- 

served chains on Si(553)-Au surface have the periodicity 
equal to the Si lattice constant in [110] direction. 

The photoemission spectrum of the Si(553)-Au surface 
is dominated by two one-dimensional bands (Si and S2) 
with parabolic dispersions [4|, [5|, [fl, LL3[ . Those bands show 
different dispersions, and as a result cross the Fermi en- 
ergy E F at different fey, i.e. at 1.04 (1.07) A -1 , and at 
1.18-1.22 (1.25-1.30) A" 1 0,0, Eft The less dispersive 
S2 band is split by 85 meV, and the splitting has been 
identified to be a spin splitting induced by the Rashba 
spin-orbit (SO) interaction [13[. Similar doublet of the 
proximal bands is also observed in Si(557)-Au surface 
[27}, which is also spin split due to SO interaction [28( | . 
The origin of the SO split S2 band comes from the hy- 
bridization of the row of Au with neighboring Si atoms 
on Si(557) terrace. Therefore it is also expected that, in 
the case of Si(553)-Au surface, situation will be similar. 
In fact the spin-split S2 band has been identified by DFT 
calculations in one of the models of Ref. [20[. However 
other features of calculated band structure for this model 
disagree with experiment. The other, more dispersive, 
band (Si) does not suffer from the SO interaction and 
is very similar to the band observed in Si(335)-Au sur- 
face Q, which also comes from the hybridization of the 
Au row with neighboring Si atoms [29[. The fact, that 
in case of Si(335)-Au surface one does not observe SO 
splitting (or SO is very small) may come from different 
widths of terraces, and thus interaction of Au row with 
the step edge Si atoms. 

Early experimental and theoretical investigations as- 
sumed Au coverage to be 0.24 ML, which is suffi- 
cient to form a single row of Au atoms per terrace 

H S S 0, H S M m m (3 Q M dzl EI . Thus 

it is no wonder that DFT calculations did not reproduce 
photoemission spectra. Recently determined Au cover- 
age on Si(553) surface is 0.48 ML, i.e. two Au chains per 
terrace [30j. In fact two Au chains per terrace are con- 
sistent with the experimentally obtained coverage from 
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x-ray diffraction [15[ . However calculated band structure 
for the structural model deduced from this experiment 
does not agree with the measured photoemission spectra 
[H, [lj[. The widely used coverage of single Au chain 
for Si(553) surface traces back to the initial assumption 
of 2/3 ML coverage for Si(lll)\/3 x V^-Au and to 0.44 
ML for Si(lll)5x2-Au reconstructions jsilsl. Thus in 
light of this new value of Au coverage, most of structural 
models of Si(553)-Au surface need to be revised in order 
to take into account two gold chains per terrace. 

The purpose of the present work is to determine a 
structural model of the Si(553)-Au surface, which accom- 
modates two gold chains, and calculate corresponding 
band structure. The structural model derived from first 
principles density functional calculations features two Au 
rows running parallel to the step edges and located in the 
middle of terraces. The step edge Si atoms rebond in or- 
der to form honey-comb structure, which is also present 
in other vicinal Si surfaces [4|, [29[ . This structural model 
is more stable than other structural models with sin- 
gle Au cha in OH [20[ and the model deduced from x-ray 
diffraction [l5[. Moreover, the calculated band structure 
for this model perfectly matches the measured ARPES 
spectra, showing two metallic bands associated with the 
hybridizing Au rows with the neighboring Si atoms. The 
rest of the paper is organized as follows. In Sec. [Ill the de- 
tails of calculations are provided. The structural model of 
Si(553)-Au surface is presented and discussed in Sec. IIII[ 
while Sees. HVl and IVl are devoted to the simulated STM 
topography images and the electronic band structure, re- 
spectively. Finally, Sec. IVII contains some conclusions. 



II. DETAILS OF CALCULATIONS 

The calculations have been performed using stan- 
dard pseudopotential density functional theory and linear 
combination of numerical atomic orbitals as a basis set, 
as implemented in the SIESTA code [H, H H, |37j . 
The local density approximation (LDA) to DFT [381 ], 
and Troullier-Martins norm-conserving pseudopotentials 
[39[ have been used. In the case of Au pseudopotential, 
the semicore bd states were included. A double-^ polar- 
ized (DZP) basis set was used for all the atomic species 
[34l [3^ . The radii of the orbitals for different species 
were following (in a.u.): Au - 7.20 (5d), 6.50 (6s) and 
5.85 (6p), Si - 7.96 (3s), 7.98 (3p) and 4.49 (3d), and H - 
7.55 (Is) and 2.94 (2p). A Brillouin zone sampling of 12 
nonequivalent k points, and a real-space grid equivalent 
to a plane- wave cutoff 100 Ry have been employed. 

The Si(553)-Au system has been modeled by four sili- 
con double layers and a vacuum region of 19 A. All the 
atomic positions were relaxed until the maximum force 
in any direction was less than 0.04 eV/A, except the bot- 
tom layer. The Si atoms in the bottom layer were fixed 
at their bulk ideal positions and saturated with hydro- 
gen. To avoid artificial stresses, the lattice constant of Si 
was fixed at the calculated value, 5.39 A. 



III. STRUCTURAL MODELS 

The total energy calculations of the Si(553)-Au surface 
show, like in case of other Si surfaces 0, 0, , that it 
is energetically favorable for the Au atoms to substitute 
into the top Si layer. The surface energy gain per unit 
cell is more than 1 eV, as compared to the adsorption 
of Au. Therefore, in the following, I will focus on the 
structural models of the Si(553)-Au surface featuring the 
top Si layer atoms substituted by the gold. 



A. Single unit cell 

Since the STM measurements 0, i, 0, 0, H i, ED, 
[ill [l2| show that the topography modulation along the 
chains is equal to the Si lattice constant in direction [110], 
it is natural to take for calculations a single unit cell in 
this direction. This can be also supported by the fact 
that the buckling is not observed in Si(553)-Au surface, in 
contrast to Si(557)-Au and Si(335)-Au surfaces [H,[25[. 

Almost all the pro posed structural models of Si(553)- 
Au surface 0,0, [13] feature single Au chain per terrace, 
as the determined Au coverage was twice as small as the 
actual one [30[. The only model, accounting for two Au 
chains, was deduced from x-ray diffraction experiment 
fl5lh and further investigated by DFT calculations [19|, 
[201 ] . According to that model, the Au atoms adsorbed 
near the Si step edges. However, none of the structural 
models managed to reproduce the phootemission spectra. 
Moreover, systematic investigations of structural models 
of the Si(553)-Au surface with proper Au coverage (2 Au 
atoms per Si(553) surface unit cell) show that the more 
stable models feature the Au atoms substituted for the 
top Si layer atoms in the middle of terrace. 

Figure [1] shows the most stable structural model (sin- 
gle cell model) of the Si(553)-Au surface, where the Si 
surface atoms (Sii-Siy) are labeled by numbers 1-6 and 
two gold atoms by Aui and Au2- This model is ener- 
getically more favorable than the most stable Si(553)-Au 
model with single Au row 0, [24[ , and the calculated sur- 
face energy is 16.45 meV/A 2 lower in the Au-rich limit 
[4Q( | . Some of the other models, in which Au atoms oc- 
cupy various top layer Si positions, from Sii to Sir (see 
Fig. [1] for labeling), have slightly higher energies. The 
differences are usually in the range of a few me v/A 2 , 
and the next 'best' structural model, in which the Au2 
occupies the Si2 position, has energy only 3.48 meV/A 2 
higher. The surface energies of the most stable structural 
models referred to the most stable model with single Au 
row, are summarized in Table H 

Although the energy differences between various mod- 
els listed in Table [J are rather small, there are other argu- 
ments supporting the model shown in Fig. [lj First one 
is that the model has the lowest surface energy. Second 
one concerns the honey-comb (HC) chain. It is widely 
accepted that the structural models of Au decorated vici- 
nal Si surfaces should posses the honey-comb chain at the 
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FIG. 1: (Color online) Structural model (single cell model) 
of the Si(553)-Au surface. Top panel shows side view of the 
structure, and bottom panel shows top view with marked sur- 
face unit cell. Labels 1-7 stand for silicon surface atoms (Sii- 
Si7), while Am and Au2 denote gold atoms. The dashed lines 
in bottom panel indicate the step edges. 



TABLE I: The relative surface energies of most stable struc- 
tural models of Si(553)-Au structure. The energies are re- 
ferred to the most stable model with single Au row. 



position of Au 


surface energy (meV/A 2 ) 


Am, Au2 


-16.45 


Sii, Au2 


-12.24 


Si2, Au2 


-12.97 


Si4, Au2 


-10.35 


Si7, Au2 


-9.28 



taken for calculations is equal to the Si lattice constant in 
this direction, all the models discussed above do not take 
into account the possibility of the buckling at step edge, 
by definition. This is consistent with the experimental 
observations @, i, S, 0, @, !, ED, El El- However, it is 
interesting to check if the doubling of the unit cell in [110] 
direction will leave the present model of the Si(553)-Au 
surface unchanged. 



B. Doubling the unit cell 

Doubling the unit cell along the step edges, i.e. in 
[110] direction, leaves almost all the structural models of 
the previous subsection unchanged. Moreover, it turns 
out that the models listed in Table [J have again the low- 
est surface energies. It seems that the Si(553)-Au sur- 
face prefers the symmetrical arrangement of Au atoms in 
neighboring surface unit cells. The models with asym- 
metrically arranged Au atoms have much higher surface 
energies. What is interesting, none of the most stable 
structural models leads to the buckling of the step edge 
Si atoms, in full agreement with the experimental results 

H H B 0, H B EE El El- Thus one can assume that 
the single cell model (Fig. [1]), i.e. the model with gold 
atoms in positions Aui and Au2, is a good candidate for 
structural model of the Si(553)-Au surface. However, a 
more detailed inspection shows that this is not the case. 
It turns out, that there is a model (double cell model), 
very similar that shown in Fig. [TJ in which gold as well 
as Sis atoms slightly change their positions, leading to 
dimerization of Au atoms in the rows. The correspond- 
ing model is shown in Fig. [2j The dimerization of the 




step edges. The main feature of HC structure is the pres- 
ence of a true double bond between Si2 and Si3 atoms (see 
Fig . [1]), which is responsible for stability of the HC chain 
[411 ]. However, none of the models, presented in Table U 
features the HC chain at the step edge, except the present 
best stable model, shown in Fig. [TJ Third argument, 
concerning the band structure, also supports the present 
most stable model. The calculated band structure re- 
produces the photoemission spectra of Refs. 0, H, S EH 
reasonably well. The other models, which have slightly 
relative higher surface energies disagree with the ARPES 
spectra. In particular they do not give correct behavior 
of bands near the Fermi energy. This will be discussed 
in Sec. [3 

It is known that other Au decorated vicinal Si surfaces, 
like Si(557)-Au or Si(335)-Au, show the buckling of the 
step edge Si atoms [24], [25j . This manifests itself in STM 
topography images as the periodicity doubling along the 
chains. Since the periodicity of the slab in [110] direction 




FIG. 2: (Color online) The most stable structural model (dou- 
ble cell model) of the Si(553)-Au surface. Top (bottom) panel 
shows side (top) view of the structure. Again, labels 1-7 stand 
for silicon surface atoms (Sii-Siz), while Am and Au2 denote 
gold atoms. 

Au atoms further lowers the surface energy. The en- 
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ergy gain is 13.73 meV/A 2 , as compared to the single 
cell model. It is worthwhile to note that the local ar- 
rangement of the Au and Sis atoms is the same as in 
recently proposed models of the Si(lll)5x2-Au recon- 
struction [43, H, EHf. Moreover, the dimerization 
of the Au atoms leaves the structure at the step edges 
unchanged, i.e. there is no buckling at the step edges, 
so the periodicity along the terraces still is equal to the 
Si lattice constant in [110] direction. However, one can 
notice a sort of horizontal buckling, i.e. the Sis atoms al- 
ternate between left and right positions in the direction 
perpendicular to the steps (see Fig. [2|), with distortion 
Ay — 0.73 A. This fact is reflected in STM topography 
and band structure, as will be discussed in text sections. 



IV. STM TOPOGRAPHY 

The STM topography data of the Si(553)-Au surface 
shows one-dimensional structures, which are interpreted 
as the step edge Si atoms 0, [HI, [H, [20|. Both struc- 
tural models discussed in Sec. IIIII support this scenario. 
To further check which one of two proposed models, i.e. 
single cell and double cell models is closer to real model 
of the Si(553)-Au surface, the STM simulations within 
the Tersoff-Hamann approach [46[ have been performed. 
The results of rilled and empty state constant current to- 
pography calculated for single and double cell models are 
shown in Fig. [3] and in Fig. [H respectively. 

Figure [3] represents simulated STM topography of 4x3 
nm 2 of the same area of the Si(553)-Au surface for sam- 
ple bias U=-1.0 V (top panel) and U=+0.5 V (bottom 
panel), obtained within the single cell model. As it was 
mentioned previously, the most visible chain structure is 
associated with the step edge Si atoms. The modula- 
tion of the topography along these chains is equal to the 
Si lattice constant in [110 ] d irection, in agreement with 
STM experiments [1, H S EL 03 • Less visible structure, 
observed at both polarizations, has been identified as due 
to Si4-Aui bonds (see Fig. [1]). Similar structure has oc- 
casionally been observed in STM experiments @, 0, [H| • 

The corresponding simulated STM image generated 
from the double cell model of the Si(553)-Au surface is 
shown in Fig. [H Top panel represents filled state (U=- 
1.0 V), while the bottom one - empty state topography 
(U=+0.5 V). Similar as in the single cell model (Fig. 
[3j), the modulation of topography along the step edge 
Si chain is equal to a^jo]- As one can notice, the less 
visible structure in filled state image (top panel) com- 
ing from the bonding of the Si4 and Aui atoms, is very 
similar to that obtained in the single cell model (see top 
panel of Fig. [3|). However, due to the dimerization of 
the Aui atoms, the Si4-Aui bonds are slightly rotated 
now, which is reflected in the STM image. On the other 
hand, the empty state topography (bottom panel of Fig. 
2]) is completely different from that of single cell model. 
We observe a sort of zig-zag structure now, which comes 
from the every second Sis atom (closer to the step edge 




FIG. 3: (Color online) STM simulations of 4x3 nm 2 area of 
the Si(553)-Au surface, calculated for the single cell model. 
Top (bottom) panel represents filled (empty) state topogra- 
phy, respectively. 



on the same terrace), Aui, Si4 and Si3 atoms (see Fig. [2]). 
Moreover, the bonded Sis and Aui atoms appear higher 
in STM topography image. This could correspond to the 
experimentally observed similar structure pj ul I9L fill, \v\ . 

As we could see, both structural models, giving differ- 
ent topography images, seem to be consistent with the 
experimental results. So a question arises, how it is possi- 
ble that we experimentally observe different topography 
images, once consistent with the single cell model, and 
another time, consistent with the double cell model. To 
be more precise, the STM topography recorded at the 
same conditions shows areas of different topographies. 
Some of the surface areas feature the structure in the 
middle of terraces with a^jo] periodicity, and the other 
areas show the structures with a periodicity doubling (see 
Fig. 1 of Ref. @, Fig. 2 of Ref. Q or Fig. 1 of Ref. 
|9[). On the other hand, the chain structure associated 
with the step edge Si atoms has always periodicity of the 
Si lattice constant in direction [110]. The most plausible 
scenario for such a behavior is that both structures are re- 
alized in real Si(553)-Au surface. The structure obtained 
within the single cell calculations appears to be a high 
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FIG. 4: (Color online) STM simulations of 4x3 nm 2 area 
of filled (top) and empty state (bottom) topography of the 
Si(553)-Au surface, calculated for the double cell model. 



temperature phase, while the double cell structure, with 
lower energy, is a low temperature phase. At intermedi- 
ate temperatures, both structures can be locally realized, 
as it is evident from the experimental data [1,0, El]. More- 
over, the presence of defects, which could further stabilize 
one of the phases at intermediate temperatures, cannot 
be omitted. Nevertheless, at very low temperature, only 
the double cell structure is expected to be observable. 



V. BAND STRUCTURE 

The experimentally measured electronic band struc- 
ture features two one-dimensional bands (Si and S2) with 
parabolic dispersions 0, 0, 03 • Those bands cross 
the Fermi energy Ep at different fey, i.e. at 1.07 A -1 
(Si band), and around 1.27 A -1 (S2 band) [7]. The S2 
band, crossing the Ep near the Si(553) surface Brillouin 
zone boundary (fey = 1.27 A -1 is split by 85 meV, and 
the splitting has its origin in the Rashba spin-orbit in- 
teraction (131 ] . Similar doublet of the proximal bands is 
also observed in Si(557)-Au surface [23], which is also 
spin split due to SO interaction [28j]. It is expected that, 



also in the case of the Si(553)-Au surface, situation will 
be similar. Indeed, the DFT calculations for one of the 
models of Ref. [20[, show that S2 band is spin-split. How- 
ever other features of calculated band structure for that 
model disagree with the experimental data. Therefore it 
is natural to check the band structure calculated for the 
present models against the ARPES data 0,11,0, El. Al- 
though the SO interaction was not included in the present 
calculations, so one cannot expect to get the doublet of 
proximal S2 bands, the other features of the band struc- 
ture should be reproduced very well, unless none of the 
proposed models is correct. 

Figure [5] shows the electronic band structure of the 
single cell model, calculated in the direction [110], i.e. 
parallel to the steps. A few surface bands are marked 
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FIG. 5: (Color online) Band structure of the single cell model, 
calculated in T-K-M' direction of two-dimensional Brillouin 
zone. The T-K-M' direction is parallel to the steps of the 
Si(553)-Au surface. The atomic character of the different 
bands is indicated by using different symbols. Open (filled) 
squares mark the band coming from the hybridization of the 
Aui and Sis (Au2 and Sie) atoms (see Fig. [1]), while the open 
circles stand for the step edge Si atoms band. The energies 
are measured from the Fermi energy (Ef = 0). 

with different symbols, reflecting character of the bands. 
The band marked with open (filled) squares comes from 
hybridization of the Aui and Sis (Au2 and Sie) atoms, 
while the band marked with open circles comes form the 
step edge Si atoms. The open and the filled square bands 
are those bands observed in photoemission experiments 
[1, H, 0]- Both bands are metallic and cross the Fermi 
energy at ky = 1.00 A" 1 (Si) and 1.18 A" 1 (S 2 ), and 
slightly deviate from the experimentally determined val- 
ues: 1.04 A" 1 and 1.18-1.22 A" 1 [5] or 1.07 A" 1 and 
1.25-1.30 A" 1 0. Moreover, bottom of the Si band 
(open squares) appears to be shifted towards Ep by 0.4 
eV with respect to experimentally determined value of 
-0.64 eV [4, 5]. The band marked with open circles, com- 
ing from the step edge Si atoms is not observed in pho- 
toemission, probably due to matrix elements. In case of 
other vicinal Si surfaces, like Si(557)-Au or Si(335)-Au, 
this band is half-occupied [24], [25j]. This band is split 
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due to the buckling of the step edges, and as a result we 
have two bands: one fully occupied, and the other one 
completely empty. Here, in the case of the Si(553)-Au 
surface, situation is different, as there is no buckling now 
and the band is not split, however due to the presence of 
additional row of Au atoms, this band becomes almost 
fully occupied. 

Slightly different band structure one obtains when cal- 
culating it for the double cell model, as shown in Fig. 
[6l Again, different symbols stand for different surface 




FIG. 6: (Color online) Band structure of the double cell 
model, calculated in T-K-M' direction of two-dimensional 
Brillouin zone. Open (filled) squares mark the band com- 
ing from the hybridization of the Am and Sis (Au2 and Sie) 
atoms (see Fig. [2}, while the open and filled circles stand for 
the step edge Si atoms bands. 

bands. Similar as previously, the band marked with open 
(filled) squares comes from hybridization of the Aui and 
Sis (Au2 and Sie) atoms, while the bands marked with 
open and filled circles come form the step edge Si atoms. 
The open and the filled square bands are those bands ob- 
served in photoemission @, H, 0] • They have similar dis- 
persions as in the case of the single cell model, however 
the bands cross the Fermi energy at different kn points, 
namely at 1.03 A -1 (Si) and 1.22 A -1 (S 2 ), which are 
closer to the experimental values, at least those of Ref. 
[5[. What is also important, energy of the bottom of 
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Si band coincides now with experimentally determined 
value, i.e. -0.64 eV 0, The band structure coming 
from the step edge Si atoms is now more complicated. 
Namely, the open circle band of the single cell model is 
split now. This splitting has its origin in the dimerization 
of the Au atoms and lateral buckling, discussed in Sec. 
IIIII The band marked with filled (open) circles comes 
from the step edge Si atoms which are located closer (far- 
ther) to the Sis atoms on neighboring terrace. Similar as 
in case of the single cell model, both bands are almost 
fully occupied owing to the presence of the Au atoms. 

All the above results show that the double cell model 
is best candidate for a structural model of the Si(553)- 
Au surface, as this model is consistent with the STM 
and ARPES data and features the lowest surface energy. 
However, it is also possible that the single cell model is 
locally realized on the Si(553)-Au surface. 



VI. CONCLUSIONS 

In conclusion, new structural model of the Au induced 
Si(553) surface has been proposed. The model accounts 
for experimentally found value of the Au coverage, i.e. 
0.48 monolayer, which gives two gold chains per Si(553) 
terrace. Like the structural models of the other vicinal 
Si surfaces, the present model features the honey-comb 
chain. However, there is no buckling at the step edge. 
The simulated STM topography images show chain struc- 
ture, associated with the step edge Si atoms. The calcu- 
lated band structure shows two metallic one-dimensional 
bands, coming from the hybridization of Au atoms with 
neighboring Si atoms. The experimentally determined 
band structure as well as the STM topography are well 
reproduced within the present model. 
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